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We calculate the strange star properties in the framework of the Field Correlator 

cn : , 

^ , Method. We find that for the values of the gluon condensate Gi = 0.006 GeV and 

^' 4 

; G 2 = 0.0068 GeV , which give a critical temperature T c ~ 170 MeV at \i c = 0, 

the sequences of strange stars are compatible with some of the semi-empirical mass- 
radius relations and data obtained from astrophysical observations. 
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In the last decades, a great effort has been made to understand the properties of nuclear 

O 

matter at densities higher than nuclear densities. From heavy ion collisions and astrophysical 
observations of compact objects, many attempts have been made to determine the equation 

^ | of state (EOS) for dense nuclear matter in both hadronic and quark phases. In relativistic 

H ■ 

heavy ions experiments (RHIC and LHC in the near future), extreme conditions of pressure 
and/or temperature are created in the interface of the colliding nuclei that simulate those 
existing in the interior of compact stellar objects or in the beginning of the universe. In 
both terrestrial experiments and astrophysical observations, one of the main goals is the 
determination of the microscopic description of the EOS of dense nuclear matter in the 
framework of the fundamental theory of strong interactions (QCD). 

The current treatments have led to theoretical results which are unable to fully explain 
the observed phenomena of compact astrophysical objects. Among these treatments, the 
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Nambu-Jona-Lasinio (NJL) model and MIT Bag model [3[ have been used in the 

study of compact stars to describe quark matter at nonzero temperature and density jj- 
6[. However, in both the NJL and MIT Bag models, quarks enter in the respective EOS 
as a free quark gas with Fermi-Dirac distribution. In the EOS of the NJL model, the 
confinement is not explicitly included, whereas in the MIT Bag model the confinement is 
represented by a bag enclosing the free quarks with a constant B which gives the energy 
density difference between the confined and deconfined phases. At low temperatures, as in 
compact star interiors, the NJL and MIT Bag EOS's are used in the large density region 
where the approximation of free quarks can be considered. This is not so at all densities 
and (low) temperatures of quark matter. Quarks strongly interact submitted to a potential 
that is large when the qq distances are large. In this case, a nonperturbative method must 
be considered. 

Recently the EOS of quark-gluon plasma was derived in the framework of the Field Cor- 
relator Method (FCM) [7| . The FCM (for a review see 8]) is a nonperturbative approach 
which naturally includes from first principles the dynamics of confinement in terms of color 
electric and color magnetic correlators. The parameters of the model are the gluon conden- 
sate G2 and the qq interaction potential V\ which govern the behavior of the EOS, at fixed 
quark masses and temperature. The FCM has been used to describe the deconfinement phase 
transition (912] and an important feature of the model is that it covers the entire phase di- 
agram plane, from the large temperature and small density region to the small temperature 
and large density region. In the connection between the FCM and lattice simulations, the 
critical temperature at /i c = turns out to be T c ~ 170 MeV for G 2 = 0.00682 GeV 4 Q. 

Very recently, the app 



considered in Ref. 



ication of the FCM to the study of neutron stars (NS's) has been 
14( where the microscopic EOS provided by the FCM is assumed for 
the internal NS cores. A farther interesting issue is the application of the FCM to the study 
of the strange quark matter (SQM) made of uds quarks and a class of compact stars called 
strange stars (SS's). The existence of SQM was first conjectured by Bodmer jl5| and Witten 
. Such a strange matter at zero temperature and pressure would have a ground state 

1^ | . The possibility 



17, 



18|. In 



lower than that of the 56 F e nucleus and be possibly found in NS's [15], 
of the SS's existence has been modeled in the framework of the MIT Bag model 
contrast to this model, the FCM naturally includes the confinement which manifests itself 
in the gluon condensate and the qq interaction potential. 
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Observationally, there is growing evidences of existing stars made of SQM. Among the 



possible SS's candidates is the X-ray pulsar Her X-l discovered in 1972 
isolated X-ray source RX J1856. 35-3754 which has_ 
sources SAX J1808.4-3658 and XTE J1739-285 



19j, the nearly 



Deen studied in [201 ] . the transient X-ray 
and the atoll source 4U 1728-34. The 



general trend is that for these sources the SS hypothesis has produced mass-radius relations 



(MRR's) more consistent with the semi-empirical MRR's obtained from observations [22l. |23 1 
than those obtained from the NS hypothesis. Very recently, the mass of the millisecond 
pulsar J1614-2230 has been measured using Shapiro delay technique [24j. By comparing the 
results among different EOS, the authors conclude that the high pulsar mass can only be 
interpreted in terms of a strongly interacting quark matter. 

In the present paper we made an application of the FCM to investigate the properties of 
the SS's in terms of the main parameters of the model. Our aim is to understand the role 
of the nonperturbative dynamics of confinement in the description of the SQM with charge 
neutrality and /3-equilibrium and the effects on the MRR's of SS's. 

This paper is organized as follows. In Sec. HIl we present the basic formalism of the 
FCM for the calculations of the important quantities. In Sec. Ill II we show some properties 
of the strange matter in /3-equilibrium and charge neutrality and calculate the strange star 
properties in terms of the model parameters. In Sec. IIV} we summarize the main conclusions. 



II. BASIC EQUATIONS 



Let us summarize the main aspects of the FCM. The gluon field is assumed to be split 
into a nonperturbative background field and a (valence) perturbative quantum field a M , 



= + dp 



(1) 



According to this description, the partition function of a system of quarks and gluons is 
given by 



rSLA 



N / £>0exp 



0=1/T 



(2) 



dr I d 3 x L tot (x, r ^ 

where denotes all fields a M , and L tat is the QCD Lagrangian density [8], and N 

is a normalization constant. In Eq.(J2]), < ••• >b means averaging over nonperturbative 
background fields B^. 
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In the FCM approach, the confined-deconfined phase transition is dominated by the non- 
perturbative correlators. The dynamic of deconfmement is described by Gaussian (quadratic 
in F^ U F^ U ) colorelectric and colormagnetic gauge invariant Fields Correlators D E (x), Df(x), 
D H (x), and D± (x). The main quantity which governs the nonperturbative dynamics of de- 
confinement is given by the two point functions (after a decomposition is made) 

g 2 (ir f [E i (x)^(x,y)E k (x)^(y,x)]) = 5 lk [D E + D? + uf 



g 2 ( tr f [Hi(x)$(x, y)H k (x)<S>(y, x)\ 



B 



B 



5 ik [D H + Df + ui^—^] - u iUk 



,0D E dDf 

dDf_ 

du 2 



H 



(3) 
(4) 



(5) 



where u = x — y and 

$(x, y) = Pexp 

is the parallel transporter (Schwinger line) to assure gauge invariance. 

In the confined phase (below T c ), D E (x) is responsible for confinement with string ten- 
sion a E = (1/2) J D E (x)d 2 x. Above T c (deconfined phase), D E (x) vanishes while Df(x) 
remains nonzero being responsible (toghether with the magnetic part due to D H (x) and 



Di(x)) for nonperturbative dynamics of the deconfined phase. 



nonperturbative part of Df(x) is parametrized in the form |8l. l25j] 



n lattice calculations, the 



D E {x) = Df(0)e 



-MM 



(6) 



where A = 0.34 fm (full QCD) is the correlation length, with the normalization fixed at 
T = fi = 0, 



7T 



D E (0) + D?(0) = -G 2 



(7) 



where G 2 = 0.012 ± 0.006 GeV 4 is the gluon condensate [26]. Another important quantity 



considered in the FCM is the large distance static qq potential given by 



mp 

f 



25] 



i>p POO 

Vi= dr(l-rT) idiD E {y/W+^ 2 ) 
Jo Jo 



(8) 



The generalization of the FCM at finite T and /x provides expressions for the thermo- 
dynamics quantities where the leading contribution is given by the interaction of the single 

eiark and gluon lines with the vacuum (called single line approximation (SLA)). From Ref. 
] and standard thermodynamical relations {^J , we explicitly rewrite in more suitable forms 
to our calculations the equations (for one quark system, Nf — 1) of the pressure 



1 2N r 



3 (2tt) ; 



(9) 
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the energy density 
2N, 



.SLA 



d 6 k 



E ~ T (T^- + m 



dJ 



E 1 



dT ' nq dfi q 

and include the number density of the quark system 

2N r 



L SLA (T,J E ^ q ) + f" LA (T,J E ^ 



n 



SLA 



d 6 k 



(2tt) 3 
dJ E 2N, 



f q SLA (T,J E , N )-f q ^(T,J E , N ) 



SLAi 
g 



T 



dn q (2vr) 



d 3 k 



f" LA (T,J E , N ) + f q SLA (T,J E ,^ q ) 



where 



(10) 



(11) 



C A (T,» q ,r 



and f" LA {T^ qi J E ) 



(12) 



E = ^Jk 2 + m|, /3 = 1/T, and J £ = V\/2T is the Polyakov loop exponent. In order to give 
Eqs. (l§j)- f[TTj) in its most general forms, we assume that V\ is, in principle, a function of 
temperature and chemical potential. However, according to the parametrization gtven by 
Eq. (jSJ), V\ does not depend on the chemical potential. As pointed out in [9|, the expected in- 
dependence of V\ should be weak for values of fi much smaller than the scale of vacuum fields 

which is of the order of ~ 1.5 GeV) and is partially supported by the lattice simulations 

28]. As in [k], [jjjl, in this work we take V\ constant. 



The pressure and energy density of gluons are given by 



SLA 
Pgl 



(27T) 



d 3 k- 



k 



e /3{k+TJE) 



and 



.SLA 



dT E 2 



dT 



(27T) 



d 3 k- 



e P(k+TJ E ) 



(13) 



(14) 



where J E = \ J E is the Polyakov loop exponent in the adjoint representation. 



In order to take into account the presence of electrons to keep the quark matter in j3- 
equilibrium and with charge neutrality, we also include the equations for the pressure, energy 
density and number density of electrons given by 



1 2 



Pe 



3 (2tt) 3 
2 



d 3 k— [/ e (T,/x e ) + /e(T,/x e )] 



d 3 kE e [f e (T,n, e ) + f e (T,fx e )} 



(15) 
(16) 
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and 



where 



n e = j d 3 k[f e (T,fi c ) - f e (T,fx e )) , (17) 



/e(T,// e ) - e p ( E D -^) + 1 ' 7e(T,/i e ) - + ~ (18) 



and E e = y^T 

Inside a SS, the composition of quark matter at a given baryon density or chemical 
potential must be in equilibrium with respect to weak interactions and the overall charge 
neutrality must be maintained. The weak interactions reactions are given by 

d -> u + e + v e (19) 

and 

s -> u + e + u e . (20) 

The above reactions result in an energy loss by the star due to neutrino diffusion. In this 
case, neutrino chemical potential may be set equal to zero and the chemical equilibrium is 
given by 

Aid = /i« + /i e (21) 

and 

Hs = Hd- (22) 
The overall charge neutrality requires that 

±(2n s u LA -n s d LA -n s s LA )-n e = (23) 

where n« is the number density of particle i. By numerically solving Eqs. (l2"T|) - (j2"3"|) . for each 
value of the input total number density 

n = n s u LA + n s d LA + n s s LA + n e , (24) 

the unknown chemical potentials /i u , /id, /i s and fi e are determined for fixed values of T, G 2 
and V\. We here also include two useful quantities to express some of our results, namely, 
the baryon number density % = (n^ LA + n d LA + n^ LA )/3 (also used in Eq. (1301 ) and the 
baryon chemical potential hb = (a*u + + / i s)/3. 
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The total pressure and energy density of the quark-gluon system, including the electrons 
inside the star, are given by 

P = P S 9l LA + E P" LA ~ A l^-l + P* > ( 25 ) 

q=u,d,s 



and 



where 



sla + ^ e s q LA + A\e mc \ + e e , (26) 



'9l 

q=u,d,s 



11 _ 2jy 

AKacI = ^^ AG 2 , (27) 



is the vacuum energy density difference between confined and deconfined phases in terms 
of the respective difference between the values of the gluon condensate, AG 2 = GziT < 
T c ) — G 2 (^ > T c ) ~ \G2 and Nf is the number of flavors. 

We use the quark masses m u = 5 MeV, rrid = 7 MeV, m s = 150 MeV, the nuclear 
saturation density n = 0.153 fm -3 , and the corresponding nuclear energy density e$ = 
0.141 GeVfm -3 . With exception of the phase diagram calculation, our investigation in the 
present work is made for T = and V\ constant. Then, in Eqs. (TlOl) . (fTTl) and (HM we have 
T 2 dJ E /dT = -V x /2, 8J E /dfi q = and T 2 dJ E /dT = -9VJ8, respectively. 

The masses and radii of sequences of strange stars are calculated by numerical integration 



of the Tolman-Oppenheimer-Volkov 29( equilibrium equations for the mass and pressure 
given by 



and 



—M(r) = 4ir 2 c(r) (28) 
dr 

^_ M{r)e{r) [1 + p{r)/e{r)} [1 + A^p{r)/M{r)] 

dr P[T) ~ r 2 l-2GM{r)/r ' 1 J 

where e is the total energy density and G is the gravitational constant. For a given central 

value of e, the numerical integration gives the gravitational mass M and the stellar radius 

R. 

In our calculation, stability of the stars of a given sequence against disassembly of the 
equivalent number of neutrons to infinity is taken into account by comparing the gravita- 



tiona. mass of each star calculated from Eq.M with the baxyon mass M A given by 30] (see 
also (3l| for details) 

r R 

M A = 47rm n n B (r)[l - 2GM(r)/r}- 1/2 r 2 dr . (30) 
Jo 
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The stability against dispersion to infinity is always fulfilled along the sequence if M < Ma- 
Points for which M > M4 in the star sequence are excluded. 

III. RESULTS 

We concentrate our results on the main features of the application of the FCM to the 
calculation of SS properties. We show some aspects of SQM with charge neutrality and in 
/3-equilibrium and investigate the SS properties in terms of the model parameters. 

A. Electrically neutral strange matter in /3-equilibrium 

In the framework of the FCM, the properties of the electrically neutral SQM in fl- 
equilibrium are determined by the model parameters G2 and V\ . In Fig. [TJ we plot the 
energy per baryon minus the nucleon mass as function of Hb/uq for different values of G2 
and V\. We note that e q /ub — rriN becomes larger as G2 increases. We also note that 
Eq LA /riB — m N increases as V\ increases because \i q (q = u,d,s) behaves as an increasing 
function of V% at fixed T = and n s q LA . Such a feature which is expressed in Eqs. (fT0i) -( jT2|) 
can be explored by using the values of G2 and V\ to verify the Bodmer-Witten conjecture 



15l . 1 1 €31 ] . We come back to this subject in Sec. IIII CI 

Fig. [2] shows the sensitivity of the EOS to the values of G2 and/or V\ where the pressure 
is depicted as function of the energy density. In this figure, we see that the EOS becomes 
softer as G2 increases. From panel (a) to (b) we also observe a slight softening of the EOS 
when V\ becomes larger. These features are explicitly expressed in Eqs. (J9j), (|T0l) . and (|T2l) 
and represent the overall nonperturbative effects of color confining forces. At this point it 
is instructive to calculate the speed of sound v s /c = (dp/de) 1 ^ 2 . For the values of G2 and V\ 
considered in this work we have obtained v s /c <, 0.577, the upper bound corresponding to 
p = e/3. 

To proceed our investigation, some remarks concerning the calculation of SS masses must 
be made. In this work, we are considering pure quark stars made of quarks u, d, s. So, 
in order to avoid the quark-hadron phase transition when the density of the SS decreases 
from its maximum at the center of the star to its minimum at the star surface, we must 
chose the appropriate values of G2 and V\ at a given temperature. In other words, G2 
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and V\ must be constrained by density values for which quark matter is deconfined and 
by allowed astrophysical possibilities for masses and radii. To this end, we first obtain the 
phase transition curve T c [p) by using the condition p{T c ) = p#(T c ), where pu is the pressure 
of the confined (hadronic) phase. Quark stars with crust have been currently considered in 
the literature. Probably, in our calculations, a hadronic crust would affect the values of 
the parameters G2 and V\. As long as the author knows, up to date a hadronic EOS has 
not been derived in the framework of the FCM as done for the quark-gluon plasma phase 
[3]. The inclusion of an EOS obtained from a theory with diferent parameters would induce 
us to misunderstand the true role of G2 and V\. In order to avoid wrong interpretation of 
the FCM parameters, we neglect the possible existence of a hadronic crust no matter how 



thin it may be. As done in 9_|, Il0| . we work in the pn = approximation. Then, from Eq. 
fl25|) we obtain the phase diagram of the SQM with /3-equilibrium and charge neutrality, by 
numerically solving the equation 

P" LA + - A l e -I + Po = , (31) 

q=u,d,s 

together with Eqs. ( l2TT) - (|24l) . to determine \i u , fid, Us, He and T for given values of n, at 
fixed G 2 and V\. 

In Fig. [3], T is depicted for different values of G 2 and V\ as function of n B /n in panel 
(a) and as function of ^b in panel (b). The general trend is that T increases with the 
increase of G2 and/or V\ for fixed tib^e})- In panel (a), we observe that at T = the curves 
corresponding to different values of V\ converge to the same value of % /no depending only 
on G2 1 . On the other hand, it is not exactly known the density at which quark-hadron 
phase transition takes place. Due to this fact, we extrapolate our calculation to obtain 
G 2 = 0.00136 GeV 4 in order to give n B /n = 1 at T = 0. We use this value of Go 
reference one to guide us to examine the results for SS masses and radii . We anticipate 
that for G2 > 0.00136 GeV 4 the results naturally show that ub/uq > 1 inside the stars, as 
well as on its surfaces. Panel (b) shows the temperature as function of \ib (for the same 
values of G2 and V\ of panel (a)). For V\ = 0.5 GeV and G2 = 0.006 GeV 4 , the phase 
transition curve T c (fi) of the electrically neutral SQM in /3-equilibrium is nearly the same 



1 From Eq. (TT2]) we have [i q = ^Jk 2 + m\ + Vi/2 for T — > 0. Thus, in Eqs. (HJ) and (J22]), Vi/2 is canceled 
(assuming V\ independent of the chemical potentials of u,d,s quarks, of course). 
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that of jsl [lo| where T c (0) ~ 0.170 GeV for rif = 3. The role of the electrons is to maintain 
the /3-equilibrium and charge neutrality among the population of u,d,s quarks. Along the 
curve T c , the electron number density is small, varying from n e /riB ~ 10 -2 at \ib — to 
n e /n B < 10~ 4 at \l b ~ 0.25 - 0.6 GeV. 



B. Strange star calculations 

Our next step is the calculation of gravitationally stable non-rotating SS sequences for 
several values of G 2 and V\. We divide our investigation in two steps. First, independently 
of the results obtained in (9), Q by comparison with lattice predictions, we construct a 
scenario for Go and V\ to calculate the main properties of SS. The calculation must here be 

nn 

taken only as an extrapolation beyond the results of [9|, H0] to obtain values for G 2 and V\ 
in order to calculate masses and radii of stable sequences. Second, we limit ourselves to the 
values of G2 in 0, Q as shown in Fig. |9] 

We now calculate SS sequences at T = from Eqs. ( 1281) . (1291) and ( 1301 for fixed values of 
the FCM parameters. In order to yield stable sequences, we constrain the values of G 2 and 
V\ with respect to the condition M < Ma and stability against gravitational collapse to a 
black hole. The former is related to the binding energy of the stars and is responsible for 
the left end points of the sequences while the later determines the right end points at the 
maximum mass. In Fig. H]we plot the mass M (in units of the solar mass M & = 2 x 10 33 g) 
as function of the central energy density for different values of G2 and/or V%. Each sequence 
begins at a minimum mass and ends at a maximum allowed mass; beyond this mass, the 
stars are no longer stable against gravitational collapse to a black hole. Between the left 
and right end points we have M < Ma (the sequences for Ma/M q are not shown in panels 
(a) and (b)). 

In Fig. we plot the respective radii, where the right end points correspond to the 
maximum masses shown in Fig. HJ In both Figs. H]and|5l the masses and radii decrease as G2 
and/or V\ increases because of the softening of the EOS. It can also be noted that the lengths 
of the sequences are shortened, with the left end points going up to the respective right end 
peaks. We explore this behavior to determine G2 and V\ in order to obtain sequences in 
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the limit of only one star at the maximum mass 2 M = Ma, as the one represented by the 
small open square dot (where the excluded parts of the sequences for which M > Ma are 
depicted to guide the eye), shown in panel (c) of Fig. HI Then, for G 2 = 0.006 GeV 4 and 
V\ = 0.0671 GeV we obtain a lower bound for the maximum masses of the sequences obeying 
the M < Ma condition. For V\ > 0.0671GeV, we have sequences with M > Ma (which 
have been discarded). We remark once again that the maximum masses diminish with the 
increase of G<i and/or Vi, satisfying the M < Ma requirement up to the lower bound given 
by M = Ma- In doing this for several values of Gi and V\, we obtain a curve in the plane 
of the FCM parameters where each point (G2, V\) represents a SS sequence in the limit of a 
single star at M — Ma-, as shown in Fig. [6j Below the M = Ma curve, all masses in the SS 
sequences satisfy M < Ma', otherwise, all sequences correspond to M > Ma- 

An important point to be observed along the M = Ma curve is that the baryon number 
density on the SS surface increase with the increasing value of G<i- In particular, we have 
found that ub/uq = 1 for G2 = 0.00136 GeV 4 . This value of G2 is too low with respect to 
lattice predictions to be realistic. Then, greater values of the gluon condensate favor SS's 
surface densities higher than the nuclear ones. On the other hand, the values of V\ become 
lower along the M = Ma curve. Also shown in Fig. [6] is the point at G2 = 0.00682 GeV 4 ( 
AG 2 = 0.00341 GeV 4 in Refs. for which SS sequences may be obtained for 0<Vi< 

0.0483 GeV. It is worthwhile emphasize that differently from the lattice prescription which 
gives V\ ~ 0.5 GeV (si, [lQ|, this value can not be achieved by the astrophysical calculations 
of the present work, except for unrealistic low values of G2 with respect to the lattice results. 
For instance, extrapolating our calculation for G2 = 0.0005 we obtain V\ ~ 0.33 GeV. We 
note that for < G2 < 0.00924 GeV 4 astrophysical calculations constrain the values of V\ 
to be < V x < Vi(M = M A ) ■ 

Fig. [7] shows the masses (in panel (a)) and radii (in panel (b)) for each point (G^, V"i) 
of Fig. [6l as well as the M < Ma and M > Ma regions. In both panels, along the 
M = Ma curves we have masses in the range 1 < M/M & < 2 (panel (a)) and corresponding 



2 The stars of an equilibrium configuration, obtained from Tolman-Oppenheimer-Volkov equations, pass 
from stability to instability (with respect to radial modes of oscillations) at a given central energy density 
£ c when the mass is stationary, e. g., dM(e c )/de c = 0. Moreover, the baryon number Na given by the 
integral in the r.h.s. of Eq. (|30| is also stationary at the same value of e c as the mass M(e c ) [3^]. Then, 
it follows that M(e c ) and Ma{£ c ) — rn n NA(s c ) are maxima at the same energy density (for a detailed 
analysis of stability, see 
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radii in the range 6 km< R <12 km (panel (b)). Additionally, for G 2 = 0.00136 GeV 4 
and V% ~ 0.24 GeV we have M/M & = 1.76 and R = 11.2 km. In the other extreme, 
for G 2 = 0.00927 GeV 4 and V x = 0, we have M/M = 1.11 and R = 6.1 km. Finally, 
for G 2 = 0.00682 GeV 4 we have M/M Q = 1.2 and R = 6.7 km, showing that for the 

nn 

prescriptions of G 2 given by [9J, H0[ astrophysical masses and radii are allowed for V\ in the 
range < V x < 0.0483 GeV. 



1. Mass-radius relations 



Although it is very difficult to distinguish SS's from NS with respect to many observable 
phenomena, their MRR's present striking qualitative differences. These differences are well 
illustrated in Fig. 3 of [33] which shows MRR's for several hadronic and SS star models. 
NS generally have masses that increase from a minimum value with decreasing radii. In 
contrast, there is no minimum mass for SS's and for M < 1M , M oc R 3 . On the other 
hand, for masses in the range 1M < M < 2M the radii of SS's can be as large as 10 km. 
However, this conclusions are valid for the MIT Bag model which has been largely used. 
This is not the case when we calculate MRR's in the framework of the FCM, as we show in 
what follows. 

To show the sensitivity of the theoretical predictions with respect to the FCM parameters 
we calculate straightforwardly the MRR's relations as done for M — e and R — e relations 
in previous sections. In Fig. El we display the variation of the MRR's for different values of 
G 2 and/or V\. The region for which R < 2GM is excluded to satisfy the condition that the 
radius R must be larger than the black hole surface radius 2GM. The curves end near the 
maximum masses, where instability against gravitational collapse to black hole sets in. In 
the left panel, depending on the values of G 2 (with V\ = 0) the masses can be as large as 
2M with radii around 11 km. By switching on V\, as shown in the right panel, we observe 
a similar trend to the case V\ = 0, but with lower masses and radii. Additionally, we observe 
that in contrast with the MIT Bag model the sequences for the FCM present lower limits 
which rise when G 2 and/or V\ becomes larger. 

Before we compare the FCM theoretical predictions with observations, let us summarize 
the characteristics of some observed sources taken as SS candidates. One of the best studied 



X-ray pulsars is Her X-l discovered by Tananbaum te al. 19| in 1972. This source is 
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generally classified as a Low Mass X-ray Binary (LMXB) with M = 0.98 ± 0.12 M and an 
estimated radius R — 6.7 ± 1.2 km [341 ] ■ 

Another interesting strange star candidates is the nearby isolated X-ray source RX 
J1856. 35-3754, from the deep Chandra LETG+HRC-S observation, with inferred mass 
M = 0.9 ± O.2M and radius R = km [35]. This source has been object of contro- 
versial claims about its nature. One group 



another group 



] argue that it is a strange star whereas 
361 ] asserts that it is a normal neutron star. The controversy arises mainly 



due to the distance estimate. The main question concerns with the measurements of angular 
diameter, and thus the radiation radius. Although it may be premature to conclude about 
the nature of this source, as pointed out in 37J], we here follow the line of |20J where RX 
J1856. 35-3754 is best described as a SS. One way to obtain information about the MRR is 

n 

by the determination of the radiation radius of the source (see [20j for details). The inferred 
radiation radius -Roo is related to the mass M and radius R of the star by 

#oo = ; R (32) 
y/1 - 2GM/R V ; 

where the factor Rj \J\ — 2GM/R comes from the gravitational redshift effect on the emitted 

radiation near the star surface. According to |20j the inferred radiation radius is R^ = 

3.8 - 8.2 km for the distance to RX J1856. 5-3754, D = 111 — 170 pc. Rewriting Eq. ((321) in 

terms of M as function of R, we can obtain the relation between the mass and radius, for a 

fixed value of R^. In Fig. [9] the MRR is represented by the dotted curves denoted by R38 

and R82. 

The source SAX J1808. 4-3658, discovered with the Wide Field Camera on board the 



BeppoSAX satellite in September 1996 38], is a transient X-ray source and LMXB. From 
measured X-ray fluxes during the high- and low-states of the source, X.-D. Li et al. 22] 
have derived an upper limit for the stellar radius in SAX J1808. 4-3658, which is given by 



r 2/7 (^) 1/3 km 



(33) 



where F max and F min are the X-ray fluxes measured during high- and low-states, respectively, 
and P is pulse period. The result is depicted in FigJHJ In [22] the results indicate that SS 
models are more consistent with the above semi-empirical MRR for SAX J1808. 4-3658 than 
hadronic star models. 

The transient X-ray XTE J1739-285 was observed by NASA Rossi X-ray Timing Explorer 



satellite 



211 ] . It is the fastest-spinning celestial body yet known with a period P = 0.891 



14 



ms. Using this rotation period, Lavagetto et al. 
star radius in XTE J1739-285, given by 



derived an upper limit for the compact 



R < 9.52(— ) 



M )1/3 
Mr. 



km 



(34) 



which is also shown in Figj9j 

The atoll source 4U 1728-34 is a LMXB observed with the Rossi X-Ra y T iming Explorer 
(RXTE) which present kilohertz quasi-periodic oscillations (khz QPOs) 40|. Based on the 



works of Osherovich ant Titarchuk 



41 



421 ]. Xiang-Dong Li et al. [23[ have derived the 



following semi-empirical upper bound for 4U 1728-34 

R< 9(^) 2/3 (^) 1/3 km 
v 1.03 y K M & > 



(35) 



where a& is given by Eq. (1) in 23], and have shown that it is possibly a SS candidate. In 
Fig. El the constraint given by Eq. (4) in [23(] is not included in the MRR of 4U 1728-34. 
Finally, the binary millisecond pulsar J1614-2230, very recently observed with the Green 



24] , which is the one with the best 



Bank Telescope, show the measured mass 1.97±0.04 M 
determination obtained to date. However, it is described as a neutron star with an interior 
region, at a given baryon density, where the transition from nuclear matter to quark matter 
takes place [43], rather than properly a strange star candidate. In the framework of our 
calculations, the mass of the J1614-2230 pulsar could only be achieved for G 2 ~ 0.003 GeV 4 
as shown in Fig. [8] (curve labeled SSI), which is a very low value of G 2 with respect to that 
obtained from lattice calculations 0,(3]. Given that it is hard to describe the J1614-2230 



24j . this subject must be carefully 



pulsar by hadronic or hybrid models, as in Fig. 3 of 
investigated. However, this is not the scope of the present paper. 

We now compare the theoretical predictions for the MRR's of SS's, for several values of the 
FCM parameters, with the semi-empirical MRR's extracted from astrophysical observations. 
The main resutls of our calculation are depicted in Fig. [9] by red, b 
which show the MRR's for Go = 0.006 GeV 4 and Go = 0.00682 GeV 4 



ue and green curves, 
[lo| and for several 



r 2 = U.UUU Vore V clilU. (jr 2 

values of the potential V\. Generally speaking we can conclude that the FCM gives sequences 
of SS's within the allowed regions of the semi-empirical MRR's. In both panels, the minimum 
masses span over a range between ~ 0.6 M Q and ~ 1.12 M Q , and the maximum masses 
between ~ 1.2 M and ~ 1.36 M & . The respective radii span over ~ 6.5 km and ~ 7.8 
km. For G 2 = 0.006 GeV 4 (panel (a)), the red and blue curves cross the error bars of 
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Her X-l and RX J1865. 35-3754 indicating that sequences of SS's are possible for values 
of Vi in the region < V x < 0.05 GeV. For G 2 = 0.00682 GeV 4 (panel (b)), the region 
of possibilities comes into the more restricted region < V\ < 0.04 GeV. By increasing 
the values of G 2 , the curves approximate to the central points of the error bars but the 
sequences become more and more restricted to lower values of V\ due to the shrinking of 
the curves. The extreme limit (although unrealistic with respect to the lattice calculations) 
corresponds to G 2 = 0.009278 GeV 4 and Vi — (see Fig. [6]) for which we obtain (see 
also Fig. [7]) M — 1.11M and R = 6.1 km, just in between the upper part of the error 
bars of the sources Her X-l and RX J1865. 35-3754. This is the lower limit for maximum 
masses that can be calculated for all possible values of G 2 and/or V\ in Fig. [61 We observe 
ihat the results of the FCM are consistent with the theoretical MRR's for SS's shown in 
4, 33, 37, 44], |45| . Then, from the above, we conclude that for values of G 2 obtained by 
comparison with lattice calculations on the critical temperature in {9], [uj] it is possible to 
calculate strange star sequences within the regions allowed by the semi-empirical MRR's 
extracted from astrophysical observations and slightly consistent with the masses radii of 
Her X-l and RX J1865. 5-3754. 



C. The Bodmer-Witten conjecture 

Since long ago, the question about the existence of an absolute ground state lower then 
that of normal nuclear matter has been considered . Bodmer and Witten pointed out 
a possibility that strange matter might be absolutely stable with respect to 56 F e at zero 



temperature and pressure [15j,|l6|. A condition for the energy per nucleon is that e/riB < m N 
and a more stringent one is that strange quark matter be energetically preferred to the 56 F c 
nucleus, which is the absolute ground state of cold matter at zero temperature and pressure, 

e/n B < M( 56 F e )/56, (36) 

where M( 56 F e ) is the mass of the 56 F e nucleus. 

In the presente work, for some choices of the values of G 2 and V\, the energy per baryon 
can be lower than the mass of the nucleon (~ 939 MeV) and the energy per nucleon in the 
56 F e nucleus (~ 930 MeV). In panel (a) of Fig. [10]we plot e s q LA /n B - M( 56 F e )/56 as function 
of V\ for several values of G 2 . Due to the linearity of the energy per baryon as function of Vi, 
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we have taken only a few points for the curves with constant G2. Depending on the values 
of G 2 , two possibilities are allowed. For G 2 > 0.0041GeV 4 we have e s q LA /n B > M( 56 F c )/56 
for all Vi > 0. When G 2 < 0.0041GeV 4 , we observe that e s q LA /n B < M( 56 F c )/56 for 
< V x < V} 0) and e s q LA /n B > M( 56 F e )/56 for V x > v} 0) , where v£ 0) is the point at which 
£q /ub — M( 56 F e )/56 crosses zero. We explore this fact to investigate the possibility of 
the existence of absolutely stable quark matter with respect to 56 _F e in strange stars. 

To this end, for each point (G 2 , Vi) of the M = M4 curve of Fig. [61 we calculate the 
energy per baryon at T = p = 0, which corresponds to quark matter conditions on star 
surfaces. The results are shown in panel (b) of Fig. [TUl where e q LA /n B — M( 56 F e )/56 
is plotted as function of V±. As in the panel (a), to each point of the M = Ma curve 
corresponds a line of constant G 2 (not all shown in the right panel (b)). All strange star 
sequences for which M < Ma are inside the triangle. We see that strange quark matter 
is unbound with respect to 56 F e for G 2 > 0.0041 GeV 4 . On the other hand, we also see 
that £q LA /n B — M( 56 F e )/56 < is only possible for values of G 2 and V\ in the regions 
0.00136 GeV 4 < G 2 < 0.0041 GeV 4 and < V x < 0.138GeV, respectively. Thus, according 
to the FCM predictions for G 2 , absolutely bound strange matter with respect to 5e F e on the 
surfaces of strange stars seems to be not allowed. 



IV. FINAL REMARKS AND CONCLUSIONS 



To summarize, we have investigated the main features of SQM with charge neutrality 
and in /3-equilibrium and SS's in the framework of the FCM. In the FCM [?]], the dynamics 
of confinement naturally appears in the EOS having the gluon condensate G 2 and the large 
distance qq interaction potential V\ as the main parameters of the model, for fixed values 
of the quark masses and temperature. 

We have considered the general aspects of the SQM at T = 0. The structure of the EOS, 
provided by the FCM, is such that it becomes softer (in the sense that for a given density we 
get lower pressure) with increasing G 2 and/or V\. An interesting feature is that the energy 
per baryon is very sensitive to the values of G 2 and/or V\ allowing for the investigation of 



the Bodmer-Witten conjecture 



15 



,3. 

We calculate stable sequences of SS's for different values of G 2 and V\. The general 
feature is that the maximum masses of the sequences and the respective radii become lower, 
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and the length of the sequences becomes shorter for larger values of G2 and/or V\. Then, 
looking for appropriate values of G2 and/or Vi, we obtain for limiting sequences with only 
one star lower bounds for the maximum masses, which determine a curve in the G2 x V\ 
plane below which all sequences are gravitationally stable. Thus, independently of what the 
lattice prescriptions for the model parameters may be, we have obtained a set of values for 
Go and V\ (which, of course, includes those obtained by comparison with lattice calculations 
[9L [lo| ) which produces stable SS sequences. The parameters are constrained by astrophysical 
conditions to spun over the ranges < G 2 < 0.00927 GeV 4 and < V x < V{ nax , where V™ ax 
is an upper bound determined by the above curve in the G2 x V\ plane. This gives us a 
general connection that open the possibilities for comparison of lattice data with results 
obtained from astrophysical observations. 

In going farther, we analyze the sensitivity of the MRR's with respect to the variation of 
the model parameters. The MRR's become shorter and with lower masses and radii with 
increasing values of G 2 and/or V\ and, depending on the values of these parameters, it could 
reach upper masses around ~ 2M and the respective radii around ~ 10 km. In order 
to make a more specific connection between nonperturbative QCD and astrophysics, we 
compare the theoretical predictions of masses and radii of SS's provided by the FCM with 
those extracted from astrophysical observations of some SS candidates. For G2 = 0.006 GeV 
and G 2 = 0.00682 GeV 4 (extracted from lattice predictions for the critical temperature at 
zero chemical potential 0, [lo)]), we obtain stable sequences within the allowed regions of the 
semi-empirical MRR's of the sources SAX J1808.4-3658, 4U 1728-34 and XTE J1739-285, 
and slightly compatible with the mass and radius of the sources Her X-l, RX J1865. 5-3754. 
We remark that it seems possible to reconcile astrophysical results with the predictions for 



the gluon condensate in [9| 



Finally, we have considered the Bodmer-Witten [15j, |l6j conjecture in the framework of 
the FCM. We have calculated the energy per baryon of SQM with charge neutrality and 
/3-equilibrium at zero temperature and pressure for different values of G2 and V\. We have 
found that only for G 2 < 0.0041 GeV 4 the energy per baryon is less than the bind energy 
per nucleon of the 56 F e nucleus. Thus, for G 2 = 0.006 — 0.007 GeV 4 absolutely bound SQM 
with respect to 56 F e on the SS surfaces seems to be not possible, unless future developments 
show to be otherwise. 
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FIG. 1: Energy per baryon of the electrically neutral strange matter in /3-equilibrium (minus the 
nucleon rest mass tun), calculated for different choices of G2 and V\ at fixed temperature T = 0, as 
function of the baryon number density n# normalized with respect to the nuclear number saturation 
density uq. 
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FIG. 2: Pressure of the electrically neutral strange matter in /3-equilibrium, calculated for different 
choices G2 and V\ at fixed temperature T = 0, as function of the energy density. The equation of 
state limit of a massless ideal gas (p = e/3) is also shown. 
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FIG. 3: Phase diagram of the electrically neutral strange matter in /3-equilibrium calculated for 
different choices of G2 and V\. Panel (a): The temperature depicted as function of the baryon 
number density ns normalized with respect to the nuclear number saturation density no- Panel 
(b): The temperature depicted as function of the baryon chemical potential hb- 
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FIG. 4: The variation of M/M & calculated for different G2 and V\ at fixed temperature T = 0. 
In panels (a) and (b), we see the shortening of the sequences when G2 and/or V\ become larger. 
Panel (c): For a fixed value of the gluon condensate, say G2 = 0.006 GeV 4 , the variation of the 
shortening of the sequences is shown for increasing values of V\ up to the limit of only one star 
in the sequence for V\ = 0.0671 GeV, indicated by the small open square dot. The long-dashed 
(M/M & ) and dot-dashed (Ma/M & ) curves serve to guide the eye. 
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FIG. 5: The same as in Fig. [4j but for the radius R. In panel (c), the small open square dot 
indicates the limiting sequence of only one star for V\ = 0.0671 GeV. 
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FIG. 6: Values of Gi and V\ for the limiting case of stable sequences with a single star at the 
lower maximum mass limit given by M = Ma , as shown in panel (c) of Fig. [J] (small open square 
dot). Below the dotted curve all sequences are stable (M < Ma)- 
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FIG. 7: Panel (a): Limiting masses of strange stars corresponding to the M = Ma curve of Fig. 
[6] as function of G2. Panel (b): The same as in Panel (a), but for limiting radii. In both panels, 
the stable sequences lie above the dotted curves. 
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FIG. 8: The mass-radius relation for different values of G2 and V±. In both panels, the curves are 
labeled by SSI for d = 0.003 GeV 4 , SS2 for Gi = 0.005 GeV 4 and SS3 for G x = 0.006 GeV 4 . The 
solid straight line gives the Schwarzschild radius as function of stellar mass. 
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FIG. 9: The mass-radius relations for different values of model parameters compared with results 
extracted from observations. The crosses are the data of the sources Her X-l and RX J1865. 5-3754, 
respectively. The black solid curves give the upper limits for SAX J1808. 4-3658, 4U 1728-34 and 
XTE J1739-285. The dotted curves indicate the radiation radii of RX J1865. 5-3754, = 3.8 
km (R38) Roc = 8.2 km (R82). The solid straight line gives the Schwarzschild radius as function 
of stellar mass. Panel (a): FCM results for G 2 = 0.006GeV 4 and V\ = (red curve), V\ = 0.03 
GeV (blue curve), V\ = 0.05 GeV (green curve), and V% = 0.0671 GeV (small open square dot) 
for the case of the limiting sequence with only a single star. Panel (b): As in panel (a) but for 
G 2 = 0.00682GeV 4 corresponding to AG 2 = 0.00341GeV 4 of Ref.[7] and for V x = (red curve), 
V\ = 0.03 GeV (blue curve), V\ = 0.04 GeV (green curve), and V\ = 0.0483 GeV (small open 
square dot ) as in panel (a). For each of the sources SAX J1808.4-3658, 4U 1728-34 and RX 
J1865. 5-3754, the allowed regions for masses and radii are in between the Schwarzschild radius and 
the corresponding curve. 
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V, (GeV) V, (GeV) 

FIG. 10: Panel (a): Energy per baryon of electrically neutral strange matter in /3-equilibrium 
minus the energy per nucleon M( 56 i ? e )/56 of 56 -F e nucleus as function of V\ calculated for different 
choices of Gi- Panel (b): The same as in panel (a), but for the points of the M = Ma curve in 
Fig. [6l showing the regions M < Ma and M > Ma- 



